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Abstract 

Within the framework of finite element systems, we show how spaces of dif¬ 
ferential forms may be constructed, in such a way that they are equipped with 
commuting interpolators and contain prescribed functions, and are minimal under 
these constraints. We show how various known mixed finite element spaces fulfill 
such a design principle, including trimmed polynomial differential forms, serendip¬ 
ity elements and TNT elements. We also comment on virtual element methods 
and provide a dimension formula for minimal compatible finite element systems 
containing polynomials of a given degree on hypercubes. 


1 Introduction 

A framework of Finite Element Systems (FES) has been developed, in [SJ [SJ [T31 [H] to 
construct mixed finite elements, generalizing those of [laiiziE]. Cast in the language 
of differential forms, following it allows for polyhedral meshes and non¬ 

polynomial differential forms. 

A minimal compatible FES (mcFES) has three key properties: it is compatible, 
it contains certain prescribed functions, and it has the smallest dimension among all 
possible finite element systems with these properties. Each property relates to a prac¬ 
tical computational purpose: compatibility is used for the design of provably stable 
mixed methods, function containment is used for estimation of approximation error, 
and dimension minimality is used to maximize computational efficiency. 

In this paper, we do the following: 

• In Section [21 we recall the main concepts of FES. We also illustrate them with 
some comments on the mixed Virtual Element Method [S]. 

• In Section [31 we show how the dimension of a mcFES can be computed in terms 
of certain cohomology groups and how a mcFES that contains a given set of 
functions can be constructed, within a larger compatible FES. These results were 
announced, mostly in French and without proofs, in m- 

• In Section m we apply this analysis and construction process to show that: 
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(i) the trimmed polynomial spaces V^. A^, defined in m, form a mcFES con¬ 
taining Vr-iA^ on simplices; 

(ii) the serendipity spaces SrA^, defined in [1], form a mcFES containing Vr-k^^ 
on hypercubes; 

(iii) the TNT elements, defined in form a mcFES containing QrA^ on hy¬ 
percubes; 

(iv) the dimension of a mcFES on hypercubes that contains T^r-A^ can be given 
in closed form. 


2 Background on finite element systems 

The notion of finite element systems is presented in detail in [SHHlIIl] so we will 
only recall the main definitions and key results that are most relevant to this paper. 

A cellular complex is a collection T of cells T, where each cell is either a singleton or 
homeomorphic to the unit ball of some Euclidean space. The cells are subject to some 
gluing conditions that make them into a regular CW complex, but where topologists 
use continuous maps for these conditions, we require the maps to be at least Lipschitz 
continuous. We stress that in the collection T, cells of all dimensions are included. 
Typically, they are taken to be fiat-faced polytopes, all sitting in some given R", but 
this is not necessary for the theory to hold. See [Ml Definition 2.1] for details. 

Given a cellular complex 7^, a finite element system is defined as in m Definition 
2.2]. If T is a cell in a cellular complex T, we denote by €^'(T) the set of fc-forms on 
T with the following property: for any T' G T included in T (including T itself), the 
pullback of the form to T' is in L^(T') and has its exterior derivative in L^(T'). 

An element system on T, is a family of closed subspaces E^{T) C £^(T), one for 
each k €N and each T € T, subject to the following requirements: 

• The exterior derivative should induce maps: 

d: E^{T) ^ E'^+\T). (2.1) 

• If T' C T are two cells in T and ixT' '■ T' ^ T denotes the canonical injection, 
then pullback by itt' should induce a map: 

Ett, : E'^iT) E^{T'). (2.2) 

For instance the spaces £^(T) constitute an element system. A finite element system 
(FES) is one in which all the spaces are finite dimensional. 

We define E^{'T) as follows : 

i;'=(r) = {u€ 0 £;'=(T) : vT,r'er r ct^ut\t'= ut'}. (2.3) 

Ter 

In this definition, which can be interpreted as encoding a continuity property of differ¬ 
ential forms, ut\t' denotes the pullback of ut to T' by the inclusion map. 

Not all finite element systems yield good spaces E^{T). As in [TU Definition 2.3], 
we consider the following two conditions on an element system E on a cellular complex 

T: 
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• Extensions. For each T G T and fc G N, the restriction operator (pullback to the 
boundary) E^{T) —>■ E^{dT) is onto. The kernel of this map is denoted Eq{T). 

• Local exactness. The following sequence is exact for each T G T: 

0-^ K-^ E°(T) —^ E^{T) -F;d™'r(r)-^ 0. (2.4) 

The second arrow sends an element of R to the constant function on T taking this 
value. 

We will say that an element system admits extensions if the first condition holds, is 
locally exact if the second condition holds and is compatible if both hold. 

In [ITJ Proposition 2.6], it is shown that for hnite element systems admitting exten¬ 
sions, local exactness is equivalent to the combination of the following properties: 

• For each T G T, E^iT) contains the constant functions. 

• For each T G 7”, the following sequence (with boundary condition) is exact: 


0-s- E°(T) -^ F;i(r)---- E‘^'’^'^(T) -^ R-^ 0. (2.5) 


The second to last arrow is integration. 

We take it for granted that a good finite element method consists of defining a 
compatible finite element system that contains certain prescribed functions and with 
which one can somehow compute. That we want a compatible finite element system is 
justified for instance in m Section 2.4], where the underlying conditions are related 
to the existence of degrees of freedom and commuting interpolators, the main tools of 
analysis of mixed methods. Typically, one would also want the finite element system 
to contain polynomials of a certain degree, to ensure corresponding best approximation 
properties in Sobolev spaces. We also remark that in some cases it is more desirable to 
contain certain exponentials [TTl [12] . 

As was already the case in [5] , one of the basic tools of our constructions is harmonic 
extension: 

Proposition 2.1. Suppose E is a FES where each E^{T) is equipped with a scalar 
product, denoted a. We suppose that T is a cell such that i2. 51) is exact. For each a G R 
there is a unique u in £'‘^‘™^(T) such that: 


u = a and Vv G El 


[dim T—1 


(T) a{u, du) = 0. 


( 2 . 6 ) 


Fix k < dimT. Any u G E^{dT) that has an extension in E^{T), has a unique extension 
in E^(T) such that: 

Vu G £’q(T) a(dM, du) = 0 and G Eq~^(T) a{u,dv)=0. (2.7) 

An element u of E^{T) such that (12.711 holds will be called E-harmonic. The propo¬ 
sition above asserts that elements having an extension have a unique A-harmonic ex¬ 
tension. Orthogonality with respect to a will be denoted T, so that (1^ can also be 
written: 

duEdE^{T) and uEdE’^-^{T). (2.8) 

We are also interested in dimension counts, so we recall [H Proposition 2.1]: 
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(2.9) 


Proposition 2.2. Let E be a FES on a cellular complex T- Then: 

• We have: 

dimE'^iT) < diuYE^{T). 
rer 

• Equality holds in lO) if and only if E admits extensions for k-forms on each 

Ter. 

Finally we recall some results on tensor products, following [S1II3]. 

Suppose lA and V are cellular complexes, equipped with FES systems B and C, 
respectively. Consider the product cellular complex: 

r={UxV : [/€ and F e V}. (2.10) 

For U GIA and V GV, define spaces A’^{U x V) by: 

A'^{UxV) = {p*uuAp\,v : uGB‘{U), vGC’^-\V), 0 < Z < fc}, (2.11) 

where pu : U x V ^ U and pv ■ U x V ^ V are the canonical projections. 

One checks that ^ is a FES on T, and we call it the tensor product of B and C. 
This name is motivated by the fact that the formula: 

u®v=p\juApyV (2-12) 

indeed defines a tensor product. Then (12.111) can be written in compact form: 

A’{U xV) = B'{U)®C'{V). (2.13) 

One recognizes a tensor product of graded spaces. 

The important result is: 

Proposition 2.3. The tensor product construction satisfies: 

• If B and C admit extensions, then so does A. 

• If B and C are locally exact, then so is A. 

• If B and C are compatible, then so is A. 

These facts were proved in and with some simplifications in m- It is possible to 
simplify the proof of the first bullet point further by using the only if part of Proposition 

O 

2.1 Comments on virtual element methods 

The Virtual Element Method is presented for scalar function spaces in [1] . It has recently 
been extended to the mixed setting, producing H(curl) and H(div) conforming spaces 
[5]. As we see it, the mixed virtual element method does essentially two things: it 
defines a finite element system and it provides a way of computing with it, that avoids 
reconstructing all the basis functions from the degrees of freedom. As an illustration 
of the methods of FES, we provide some details on the first point. The second point 
however, we leave open. 

The constructions of [5] were based on a notion of harmonic extension for differential 
forms (both continuous and discrete). Canonical degrees of freedom for compatible FES 
can be defined, as in projection based interpolation, by non-homogeneous harmonic 
extensions - see [SI Proposition 3.23] and [131 Proposition 5.44]. As we shall see, the 
mixed VEM spaces are defined by a similar technique. 
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Lemma 2.4. Let T be a cell of dimension n. 


• Fix k < n. Choose f € L^A^(T) and g G ^{T). The system, with unknown 

uee^T): 

d*du = /, (2.14) 

d*u = g. (2.15) 


has a solution if and only if d* f = 0 and d*g = 0. In this case the solution is 
unique. 

• The case k = n. Fix f G M. and g G L^A^~^(T). The system, with unknown 
u G €’^{T): 


ju = f, (2.16) 

d*w = g, (2.17) 

has a solution if and only if d*g = 0. In this case the solution is unique. 

Proof. Left to the reader. □ 

Consider now a cellular complex T. We suppose we have, for each integer k and 
cell T G T, with k < dimT, a finite dimensional subspace Z^{T) of L^A^(T), whose 
elements g satisfy d*g = 0. For k = dimT, we put Z^{T) = 0 and note that for 
u G L^A^(T) we have d*dit = 0. We define: 

A^{T) = {uG €^{T) : V T' <1 T dMwIr' € Z’^{T') and d*u\T' G Z’^-^iT')}, (2.18) 

where T' <\T indicates that T' is a face of T. 

Proposition 2.5. The above defined spaces A^{T) constitute a compatible finite element 
system and for k < dim T ; 

dimAg(T) = dim. Z'^(T) + dim Z'^-\T), (2.19) 

whereas for k = dimT.' 

dimA^{T) = l + dimZ^-\T). (2.20) 

Proof. The definition of A^{T) ensures that it is an element system. The technique of 
harmonic extension shows that A has the extension property. The local exactness in the 
form (12.51) is also trivial to check. Thus A defines a compatible element system. The 
dimension counts follow from Lemma □ 

There are several variants of mixed VEM, but the main one seems to correspond to 
the following choice of Z. We suppose that each cell is flat, so that polynomials are well 
defined objects. Fix r > 1 and put, for k < dimT. 

Z’^{T) = {f GVr-iA'^{T) : d*f = 0}. (2.21) 

Then A, defined as above, is a compatible FES by the preceding proposition, and A^{T) 
contains polynomial /c-forms of degree r. There is a notion of degrees of freedom for 
FES. Since it is rather intuitive we omit the definition; details can be found in [H 
Section 2.4]. 
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Proposition 2.6. For any integer k and cell T of dimension n, consider the linear 
forms, defined on k-forms by: 

u / u Av, (2.22) 


for some A" ^'(T). 

These linear forms constitute unisolvent degrees of freedom on the FES A. 


Proof. Let * 
Vr-lA^-^{T) 

denote the Hodge star operator. We notice that u G Vr-i 
. Moreover, for such u we have: 

A^-(T) iff ku G 


d*u = 0 <;=> d*it = 0. 

(2.23) 

Therefore, for k < diinT = n (recalling [3J Theorem 5.4]): 



dimZ'=(r) = dim{u G Vr-iA'^~‘"{T) : du = 0}, 

(2.24) 


= dim{u G PpA”-'=(r) : du = 0}. 

(2.25) 

We also have: 

dimZ'^-^(r) = dim{u G Vr-iA'^~^+^(T) : du = 0}, 

(2.26) 


= dim{du : v G PpA"-'=(T)}. 

(2.27) 

It follows that: 

dim Z'^{T) +dim Z’^-^T) = dimPpA”-'=(T). 

(2.28) 

In the case k 

= dimT we may also check: 



1 + dimZ'=-i(T) = dimPpA”-'=(T). 

(2.29) 

Therefore: 

dimAg(r) = dimPpA"-'=(r). 

(2.30) 

Suppose now 

that u G Aq(T) satisfies: 



Vv G PpA”-'=(r) y u A u = 0. 

(2.31) 

We have d*du e P^_iA'=(r), hence d+du e Pr-iA"-'=(T) C V-A^-^iT) 

. Hence: 


0 = 


u A d+dw = ± 


dw A *dM = ± 



(2.32) 


Therefore drt = 0. 

Let f) denote the standard homotopy operator used to prove the Poincare lemma, 
as in e.g. m- Recall that f), which is proportional to k on homogeneous polynomial 
differential forms, maps PpA"“^+^(T) to PpA"“^(T). We have d*u € Vr-iA^~^{T), 
so that d*u e Pr-iA"“^+^(r) and hence f)d*M € PpA"“^'(r). We write: 


M A fid*u = — J 


0 = / M A fid*u = — / MAdf)*u+ / uA-kU, 


=±y duAf)*u+y iup=y lui^ 


(2.33) 

(2.34) 


Hence u = 0. 

It follows that the integrated wedge product is an invertible bilinear form on the 
product Aq{T) X P“A"“^(T). This concludes the proof. □ 
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3 Construction of minimal finite element systems 

The material presented in this section is an expanded version in English of material 
that appeared previously, in French and without proofs, in m- Given any sequence 

-x'=-i(r) —^ X'^{T) X^+\T) —^ • • • 

satisfying d o d = 0, we use the notation kerd|X^ to denote the kernel of the d map 
whose domain is X^ and dX^ to denote the image of that map. The fc-th cohomology 
group associated to this sequence is the quotient: 

H'=(X') = (kerd|X'=)/dX'=-\ 

3.1 Sufficient conditions for minimality 

Proposition 3.1. Suppose that A is a finite element system, and that B is a compatible 
finite element system containing A. Then we have: 

dimB^{T) > dimdg(T) + dimH'^'+i(d5(r)). (3.1) 

Proof. By the rank-nullity theorem and the definition of cohomology groups, we have 
that: 

dimdg(T) = dimddS(r) -hdimkerd|d§(T), (3.2) 

dimH'=+i(d5(T)) = dimkerd|dg+^(r) - dimddg(T). (3.3) 

Hence: 

dimA'^{T)+dimR'‘+\A'o{T)) = dimkerd|d|;+i(r) -k dimkerd|Hg(r). (3.4) 
Now, since Bq{T) is exact we may choose a subspace V of Bq(T) such that the map: 


d:H^kerd|Hg+^(r), 

(3.5) 

is an isomorphism. Then it is clear that: 


Hnkerd|Hg(r) = 0. 

(3.6) 

It follows that: 


dim (T) > dim V + dim ker d | dg (T). 

(3.7) 

Since dimP = dimkerd|dg'''^(T), we apply (13.411 to complete the proof. 

□ 

We will use this as follows: 


Corollary 3.2. Suppose that A is a finite element system, and that B is 
finite element system containing A. Suppose that: 

a compatible 

duaB^{T) = dimHg(T) -k dimH'=+i(d5(r)). 

(3.8) 


Then B is minimal among compatible finite element systems containing A. 

We will show in the next two sections that, given a finite element system A, there 
exists a compatible finite element system B containing A such that (13.81) holds. Thus 
the condition is not only sufficient, it is necessary for minimality. Moreover, in the 
case that A is included in a compatible finite element system B, it is always possible 
to find a minimal one containing A, inside B. Finally, in this context, we will actually 
construct a particular minimal FES. We will show how other authors have proposed 
different definitions of hnite element spaces, which turn out to be minimal among those 
containing certain differential forms. 
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3.2 First step 

Given a finite element system A, the first step is to construct a finite element system A 
containing A such that Aq is exact. In the case that Aq is exact already, this step can 
be skipped. 

For k < dimT choose a subspace such that: 

dAl(T)®H'^+^{T) = kerd|ylg+^(T). (3.9) 

Then choose a subspace E^(T) of Bq{T) such that d restricts to an isomorphism: 

d: E'^{T) -)■ H^+\T). (3.10) 

For k = dimT, if A^(T) contains an element with integral 1, set E^{T) = 0, if not 
choose an element in B^{T) with integral 1 and let E^{T) be its linear span. 

We remark that A^{T) n E^{T) = 0 and set: 

i'=(r) := A'=(T) 0T'=(T). (3.11) 

We remark that we have an induced isomorphism: 

d:E^{T)^B^+^{AHT)), (3.12) 

where if /c = dimT, the exterior derivative operation is replaced by integration. In fact 
we might take this, together with E^{T) C Bq{T), as the defining properties of E^{T), 
and the space H^(T) introduced above merely as a device to make it clear that such an 
E^{T) exists. 

We could also have specified E^{T) from a choice of scalar products as follows: If 
k < dimT, set: 

E’^iT) := {u G B^{T) : du G Ag+^(T), du _L dA^'(T) and u _L dT^"^(T)}. (3.13) 

If fc = dimT and if A^{T) contains an element with integral 1, set E^{T) := 0, if not 
set: 

T'=(T):={uGB'=(T) : u _L dB^"^(T)}. (3.14) 

Then one checks immediately that E^{T) is a subspace of Bq{T) such that (13.121) is an 
isomorphism. 

Proposition 3.3. Under the above hypothesis, A is a FES containing A such that the 
sequences: 

0 ^ i[](T) iJ(T) -^ i'^™^(T) K 0, (3.15) 

are exact. 

Proof. That A is stable under restrictions and the exterior derivative is immediate, so 
we only need to prove the exactness property. 

Consider first k < dimT. Suppose u G Aq(T) is such that du = 0. We have 
^o(^) = ® E^{T), so we can write u = v + w with v G Aq{T) and w G E^{T). 

Then dw = —dw G dAg(T) so dw = 0, so re = 0. 

Now we have du = 0. If fe = 0 this gives u = 0. If fc > 0, we have already remarked 
that dAg“^(T) 0 dT'^“^(T) = kerd|Ag(T), so we can write v = dt with t G Aq“^(T). 
This concludes the case k < dimT. 

If fc = dimT, the same proof goes through, replacing the exterior derivative by the 
integral for the involved fc-forms. □ 




3.3 Second step 

We now consider the situation of a finite element system A included in a compatible 
finite element system _B, such that the following sequences are exact: 

0 ^ A'^iT) Al{T) -^ ^ R 0. (3.16) 

For each cell T we will augment the spaces A’(T) in such a way that the extension 
property holds, while ^o(^) preserved (i.e. no trace-free elements will be added). We 
start by carrying out this process on ID cells, then on 2D cells, and so forth. We will use 
the following result, where tr denotes the trace (pull-back) operator B*(T) B*{dT). 

Proposition 3.4. Let I > 0 and suppose T is a cell of dimension I -|- 1. Fix k < 1. Any 
u G B^{dT) such that du G tr A^+^(T) has a unique extension u G B^(T) satisfying: 

du G A'^+^T), du ± dA^{T) and u±dB^-\T). (3.17) 

Proof. Existence. Suppose first k < 1. Let v be the A-harmonic extension of dw. Since 
dt! is A-harmonic and zero on the boundary we have dti = 0. Choose u' G B^{T) such 
that dit' = V. We have d(trM' — u) = 0. If fc = 0, trit' — m is a constant c and u' — c 
is an extension of u with exterior derivative v. If on the other hand fe > 0 choose 
u" G B^~^(T) such that dtru" = tru' — u. Then u' — du" is an extension of u with 
exterior derivative v. Adding an element of di?Q~^(T) we can ensure that it becomes 
orthogonal to di?Q“^(T). 

Suppose now k = 1. Pick v G A^^^{T) such that J v = f u and v _L dAg(T). Pick 
u' G B^ifT) such that du' = v. By Stokes we have f (tru' — u) = 0, so we can choose 
u" G B^~^(T) such that dtru" = tru' — u. Then u' — du" is, as before, an extension 
of u with exterior derivative v and by adding an element of dBg~^(T) we can ensure 
orthogonality to dBQ~^(T). 

Uniqueness. If u G Bq{T) satisfies (j3.17ll we have du = 0. The condition u T 
dilg“^(T) then gives u = 0. 

□ 

Suppose we have realized our plan for cells of dimension at most I (/ > 0), the 
augmented FES being denoted A. We consider a cell T of dimension Z -|-1. 

We equip A^{dT) with a scalar product, for instance the sum of the local ones, 
defined for T' <l dT. For k < I we put: 

F'^ = {uGA'^{dT) : uEtrA^iT) and du G tr A’^+\T)}, (3.18) 

and denote by F^ the extensions of elements of F^ defined by Proposition 13.41 above. 
We also denote: 

= {uGA'^{dT) : duTtrA'=+i(T) and uE{v GA^{dT) : du = 0}}, (3.19) 
and denote by the i?-harmonic extensions of elements of G^. Finally we put: 

ifc(T) = A'=(r)-HF'=-tG'=. (3.20) 

For fc = Z -I- 1 we take simply A^(r) = A^{T). 

Proposition 3.5. On the (Z -|- l)-skeleton ofT, A is a compatible finite element system 
intermediate between A and B such that Aq{T) = Aq{T) for all k and T. 
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Proof. We remark that dF^ C A^+^(T). If fc < /, dG^ C F^^^(T) and ii k = I, & = 0. 
Therefore dyl^(r) C We have a direct sum decomposition: 

i'=(5T) =trA'=(r)©F'=©G^ (3.21) 

It follows in particular that tr : A^{T) —>■ A^{dT) is surjective. Moreover if u G A^{T), 
V £ F^ and w £ & satisfy tr(it + u + w) = 0, then tru = 0 and trw = 0, which gives 
u = 0 and w = 0. If in fact u + v + w = 0, we will also have u = 0. Therefore the sum 
(13.201) is direct and: 

AI{T) = AI(T). (3.22) 

This completes the proof. □ 

The construction is therefore complete, and provides a minimal compatible finite 
element system containing A, by Corollarv l3.2l 


4 Minimality of some finite element systems 

We now consider more concrete examples, based on polynomial differential forms on 
simplices or hypercubes. The two steps of our general construction enter differently in 
them: 

• In Section im on trimmed differential forms, the local cohomology groups (Aj (T)) 
are non-trivial, but once this problem is fixed, the obtained spaces turn out to have 
the extension property. 

• In Section 14.21 on serendipity elements, the local cohomology groups are trivial, 
so one only needs to add spaces of forms that ensure the extension property. 

• In Section 1331 on TNT elements, the extension property holds for the spaces one 
starts with, but local sequence exactness fails. By adding tiny bubble functions 
E^{T) to fix the cohomology requirement, one also needs to add forms that extend 
them to higher-dimensional cells. 

• In Section r4.41 in which we introduce “small pleasures” elements, the sequence one 
starts with has nontrivial cohomology and does not satisfy the extension property. 

4.1 Minimality of trimmed polynomial differential forms 

We now adopt notations from finite element exterior calculus [2] . We first consider finite 
element systems on a simplicial complex T. Let VrA^iT) denote the space of A:-forms 
on T £ F whose coefficients are polynomials of degree at most r. The space V~A^{T) 
is defined with the help of the operator k on forms, which is contraction by the vector 
field X 1 -^ X (different choices of origin yield the same space). Then: 

V-A^{T) = {u£ VrA\T) : ku £ VrA^-\T)}. (4.1) 

In keeping with our previous notations, VpA q{T) is the subspace of 'P“A^(T) con¬ 
sisting of differential forms whose pullback to the boundary dT of T is 0. We also 
use freely that these spaces form exact sequences under the exterior derivative, which 
follows from the fact that they constitute a compatible FES. The underlying facts are 
proved in standard finite element language in [5] , but it is also checked with the general 
tools of hnite element systems in M- 
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Similar statements hold for the Vr^^{T) spaces, but here sequence exactness requires 
decreasing r through the complex. That is, under the exterior derivative, the sequence: 


Vr+iht\T) ^ iP,Ag(T) ^ (4.2) 

is exact, whereas the sequence: 

VrA'^-\T) VrA^{T) VrAo'^\T), (4.3) 

in general is not. The cohomology group of the latter sequence is denoted: 

U^VrAoiT). (4.4) 

Lemma 4.1. ITe have : 

dV-A^{T) = dVrA^oiT). (4.5) 

Proof. Pick u S PrAg(T) and set v := du. Then v € Vt-iAq'^'^{T) C V~Aq'^^(T) and 
dt! = 0. Hence there is w € V~Aq{T) such that dw = v. □ 


The following identity shows that the spaces A^(T) constitute a minimal com¬ 


patible FES containing the spaces T’r-iA^{T). 

Proposition 4.2. We have: 

dimP-Ag(T) = dimP^_iAg(T) -h dimH'=+ip^_iA5(T). (4.6) 

Proof. We use the vertical line notation | to denote restriction of an operator to a 
subspace. We have: 

dim'P~AQ(T) = dimd'P“AQ(r) -|- dimkerd|'P~AQ(T), (4.7) 

= dirndP^A^Xr) + dimdPpAg-^(r), (4.8) 

= dimdPrA^(r) dimdP^A|;-^(r), (4.9) 

and: 

dimPj._iAg(T) = dimdPr-iAg(T) -|- dimkerd|Pr-iAQ(r), (4-10) 

= dimdVr-iA'^iT) + dimdP^Ag-i(r). (4.11) 

Moreover: 

dimH'=+iPr_iAS(T) = dimkerd|Pr_iAg+^(T) - dirndP^-iA^jX^’): (4-12) 

= dimdPrAg(r) - dimdP^_iAg(r). (4.13) 

The result follows by addition of the above equalities. □ 

4.2 Minimality of serendipity elements 


In [T], Arnold and Awanou define spaces of serendipity hnite element differential forms 
on n-dimensional cubes, denoted 5rA^(/"). They prove a subcomplex property (The¬ 
orem 3.3) and trace property (Theorem 3.5), which is equivalent, in our notation, to 
showing that the spaces 5r-fcA^(/"’) constitute a FES, for any fixed r > n. Arnold and 
Awanou also define these spaces for r < n, but then one does not get a full complex, 
but one that stops at k = r. In many applications this is useful, since one does not 


11 


require the full complex but rather two consecutive spaces, possibly three. But the FES 
framework has not been designed for this. 

They also prove that the spaces provide unisolvent degrees of 

freedom (Theorem 3.6). By O Proposition 2.5], this guarantees that the extension 
property holds. They also check the sequence exactness property (page 1566). Thus, 
in our language, they have shown that the serendipity finite element spaces constitute 
a compatible finite element system. The inclusion T^^-fcA^)/") C is obvious 

from their construction and it is implicit from their Proposition 3.7 that Aq{I^) = 

iP,_fcA§(/"). 

Thus one can interpret many of their results, by the statement that the serendipity 
elements define a minimal compatible hnite element system containing 7^r-feA^(/"). We 
now explain this, proving the main statements above, in our framework. 

It should be noted that the construction in [Tj is very different from the one we 
propose in Section 13.31 (the first step in our construction is not necessary here, for 
reasons provided below). Our construction depends on a choice of compatible finite 
element system B containing For that purpose we could take the spaces: 

= (4.14) 

We would also need a scalar product. One can use scalar products, or more algebraic 
ones, for instance one that makes some preferred basis in B orthonormal. 

Remark 4.1. We will take I to be the interval [0,1] while Arnold and Awanou take it 
to be [—1,1]. Also, we will only address the cases r > n, although the arguments can 
be adapted to treat lower order cases as well. 

Lemma 4.3. Let B^ be a space of differential k-forms on 7”“^, for each 0 < a < r 
and for each k. Define the sum of graded tensor product spaces of forms on 

r 

A-= ^B;®iP„A-(7). (4.15) 

a—0 

Consider and the n-th trace operators, defined respectively as pullbacks by the 
injections x {0} —>■ 7” and 7"“i x {1} ^ 7". Then the intersections of the kernels 
of 6° and on A' is : 


(ker,50|A-) n (kerJ^ jA*) = ^ B'^ ® . (4.16) 

a—0 

Proof. We let Aq (resp. Ai) denote the element of ViA°{I) taking the values 1 (resp. 0) 
at 0 and 0 (resp. 1) at 1. Observe that for a > 1 we have a direct sum decomposition: 

iPaA’(7) =RAo©MAi ©iP„A5(7). (4.17) 

We also have: 

iPoA*(7) =R(Ao + Ai)+]RdAo. (4.18) 

Now take an element u of A^ and write it as a sum: 

r 

U = Vo® {Xo + Ai) + Uq ® dAo + '^{Vao © Ao + Va^ © Ai + Va), (4.19) 

Q;=l 
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with: 


(4.20) 


Vo 


G z;' G B’l-\ and ® 


Observe that if u is in the intersection of the kernels, we have: 


0 = = Vo -f ^ZZao, 

a—1 

(4.21) 

r 

0 = Slu = Vo + ^ Vai. 

a—l 

(4.22) 

r 

U = Vq 0 dAo -t- ^ Va, 

(4.23) 

a—0 



□ 


Then we are left with: 


and this completes the proof. 

Proposition 4.4. On let Xj denote the j-th coordinate function. We have: 

( U €Vr-2{n-k)iI^)y 

Pr.AQ(/”) = span < M Xj(l — Xj)da:j : JC n}, 

[ #J = k. 

where, for J = {ji,... ,jk} we define dxj = dxj^ A ... A dxj^,. 

Proof. We have: 

P,A'(/")= ^ P„,A-(/)(8)...®P„„A'(/), 

Q'i + ...+Q'7i=r 


(4.24) 


(4.25) 


where we sum over integer ai in the range [0, r]. 

In order to determine PrA5(/") we apply Lemma 14.31 or its equivalent, in each 
coordinate direction. We get: 


VrA’oin = ® ® 'Pc.Ml), 

Now, concerning each factor, we have with x the coordinate function on I: 

P,A[J(/) = {u(x)a;(l - x) : u G Pg_ 2 (/)}, 


(4.26) 


(4.27) 

(4.28) 


It follows that the left hand side of (14.241) is included in the right hand side. The 
inclusion in the other direction is trivial, so the proof is complete. □ 

Remark 4.2. The result just shown seems to be taken as self-evident in [U page 1569]. 

Proposition 4.5. The integrated wedge product defines an invertible bilinear form: 


{j -A-) ■■ VrA^n X Vr-2(n-k)A--\n 


(4.29) 
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Proof. It follows from the preceding proposition that the two spaces have the same 
dimension. 

Fix w G Vr^oil^)- We want to show that the map sending 'ci— J u A v G ^ is 
non-zero only if u is zero. Write u in the form: 


u = 


E 

j 


UjJlxjil 



with: 

UJ G 'Pr-2(n-k)[I'^), JC n} and #J = k. 

Letting J' denote the complement of J in {1,..., n} we define: 

V = ^ ejUjdxj, G Vr-2(n-k)^"-\n- 


(4.30) 


(4.31) 


(4.32) 


Then we have: 

u Av = ejMj Xj{l — Xj)dxj A dxji. (4.33) 

J j^J 

We then choose the signs ej such that ejdxj A dxj/ is the volume form on Then 
/ u A is zero iff u is zero. □ 

Proposition 4.6. Fix r > 0. 

For 0 < k < n, the following sequence, for the exterior derivative, is exact: 

Vr+iAt\n ^ VrA^n ^ P.-iAg+i(/"). (4.34) 


The sequence: 

Vr+iA^~^{r) ^ VrA'^ir) ^ R, (4.35) 

is also exact (where the last arrow is integration). 

Proof. Consider the first statement. By Proposition 0^ the sequence is paired with the 
following one: 

P,_l_ 2 („-fe)A-^-+l(J”) ^ Vr-2in-k)A--Hn ^ P.+ l-2(n-fe) A"-''”'(/") . (4.36) 

The differential in this sequence is again the exterior derivative. Recall also that: 


du Av = dz u A dv, 


(4.37) 


when (for instance) u is zero on dT. Since the latter sequence is exact, the former one 
must be too. 

Consider now the second statement. The sequence is paired with: 

P^_iA"-i(/") ^ nA°(/”) ^ R. (4.38) 

On the right we have inclusion of the constants. This sequence is exact. □ 

We may summarize our results as follows: 
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Proposition 4.7. Fix r > n. Let he a compatible finite element system on 

hypercubes of dimension at most n, sueh that: 

Vr-kA^il^) C (4.39) 

Then is a minimal compatible finite element system with this property if and 

only if: 

4(/-)=p,_fcAg(/™). (4.40) 

In this case, the spaces Vr+k- 2 mA"^~^ {I"^) provide unisolvent degrees of freedom on A. 

Proof. The first statement follows from the dimensional characterization of minimality 
provided by Corollarv l3.21 given the exactness property proved in Proposition 14.61 The 
second statement recalls Proposition 03] □ 

Again we point out that [I] and Section 0] use quite different means to define a 
compatible FES A such that (14.3911 and (14.401) hold. 

4.3 Minimality of TNT elements 

We now consider the work of Cockburn and Qiu regarding the ‘TiNiest spaces containing 
Tensor product spaces of polynomials’ on cubes or TNT elements for short [13]. We 
will show how to derive their minimality claims from the framework of FES. Moreover 
we cast their constructions in the language of differential forms, which is uniform in the 
treatment of dimensions and degree of forms. We will also make explicit use of other 
facts from homological algebra such as the Kunneth theorem. 

Fix r > 1. Define the spaces, on hypercubes: 

A'{r) = Q^A-(/") =P^A*(/)(g)...(8>P^A-(/). (4.41) 

We remark that on an interval 'PrA'{I) has the extension property. Therefore, by 
Proposition 12.31 A also has the extension property. However the sequence exactness 
fails. We want to define a minimal compatible hnite element system containing A. This 
is accomplished for dimension n = 3 in |15j . 

By 0 Lemma 3.10] we have: 

AW) = 'PrKil) 0 ... 0 VrAW- ( 4 . 42 ) 

Looking at the sequence: 

0 ^ VrA°{I) VrAl{I) 0, (4.43) 

we notice that there is a non-zero cohomology group only in the 1-forms, and that it 
has dimension 2. We write out the Kunneth theorem: 

H'AS(/") RiH-(iP,A5(/))®...®H-(nA5(/)). (4.44) 

Since, as we remarked above, H°(7^j.A5(/)) = 0, most of the terms cancel. There is only 
one non-zero cohomology group in A'{I'^), namely the one in n-forms. Explicitly: 

A^A'W) = for k<n, 

H"A5(/") A\VrAW) ® ® H^(^r.A5(/)), 
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(4.45) 

(4.46) 






and the latter space has dimension 2”. 

Now, the sequence that we want to make exact is really: 

0 ^ A°{r) Al{r) ^ A^{r) ^r^ o, (4.47) 

where the second to last arrow is integration. Thus the notation in (I4.46P is slightly 
different from the one used in the beginning of the article, e.g in (13.811 . the difference 
being whether we consider that there is an R. in the end of ^5(1") or not. We notice 
that contains elements with non-zero integral. We thus have: 

Proposition 4.8. Let B denote a compatible finite element system containing A, as 
defined in Then B is a minimal one containing A if and only if, for all k and 

n: 

4 ( 7 ") = A^oin, fork^n-1, (4.48) 

dimS”-i(/") = dimy4”-i(/") -f 2” - 1. (4.49) 

This provides the dimension count for any minimal finite element system containing 
A, using Proposition [5^ 

Any construction of a minimal compatible finite element system B containing A 
would need to add spaces £'”“^(/") to A"“^(/") such that: 

dimA”-i(/”) = 2”-l, (4.50) 

B”-1(/") = A^-\r) © £;"-!(/"), (4.51) 

as in the ‘first step’ outlined in Section 13.21 The second step, outlined in Section 13.31 
would have to ensure that the extension property holds, after these additions, since they 
also take place on the faces of 

To be more specific, we notice that if x denotes the coordinate function on the interval 
/, we can choose two generators of H^('Pj.A5(/)), to be dx and P{x)dx, where P is a 
polynomial of degree r with integral 0 on I. We also impose P(1 — x) = B(x) if r is even 
and P(1 — x) = —P(x) if r is odd (to ensure that our construction becomes invariant 
under reflections in the coordinate directions). Then we let Q be the polynomial of 
degree r + 1 such that Q{0) = <3(1) = 0 and Q' = P. 

Denote by Xi the f-th coordinate map on Let J denote a non-empty subset of 
{1,..., n}. Define the n-form fj on J": 

fj = (^P{xj)dxj(^dxj, (4.52) 

where the notation is taken to mean that we have the form P{x)dx in the directions of 
J and the form da; in the other ones (thus there is a permutation in the indices going 
on). Since J is non-empty, fj has integral 0. Notice that we have thus defined 2" — 1 
linearly independent forms in A'^{I^) with zero integral. 

We now observe that: 

fj = with gj = Q{x^) (^ P{xj)dxj (^ dxj. (4.53) 

This notation is taken to mean that we can obtain gj from fj by replacing, for each 
index i in J, the term P{xi)dxi by Q[xi), and summing over i G J. 

We define: 

P"-i(J")=span{gj : J C {1,..., n}, J ^ 0}. (4.54) 

We remark that: 
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Lemma 4.9. The elements of -^^ero on dl'^ and the map defined in \3.12\) 

is an isomorphism. 

Thus they are adequate for the first step of our general construction. Notice also 
that these spaces are invariant under permutation of indices and reflections Xi i—>• 1 — Xj. 

For the second step, one should keep in mind that we have added spaces 
for all faces T of dimension d in and that the extension property is therefore no 
longer guaranteed. We could proceed by our general method. But inspired by [15], we 
prefer the following explicit approach: 

We determine a face T of /" by a choice of indices Jt C {1,... ,n} and for the 
indices i G {1,..., n} \ Jt a number Nrii) £ {0,1}. Then the associated face T is: 

T = {(xi,... ,x„) : for i ^ Jt Xi = A^t(*)}- (4.55) 

Let T be a J dimensional face of An element of E‘^~^{T) is extended to /" as follows. 
Let u be the element in question. Its extension is chosen to be: 

u = u(g) (4.56) 

i^Jt 

We then define: 

B'=(/") = A'=(/") + ^ A'=(T), (4.57) 

dim T —A:+l 

where the sum is over fc + 1 dimensional faces of /" and E^{T) denotes the space of 
extensions of elements of E^{T) defined by (14.561) . 

Proposition 4.10. The definition i4.5T\) provides a minimal compatible FES containing 

A. 

Proof. Suppose that (14.561) holds for some u G E^{T). Then: 

du = u (g) xf"^*^(l - + 

i^Jt 

^ M (g xf^^*^(l - XiY~^'^^'^\±dXi)i=j, 

j^JT 

G A^+\I^)+ E^^\T). 

dim T—/c+2 

Thus B satisfies the requirements to be an element system. Let T be a fc + 1 dimensional 
face of /" and let T' be another one. Let’s take an element u of E^{T) and restrict u 
to T'. Choose i £ Jt\ Jt'- The restrictions of u to the faces x^ = 0 and Xi = 1 inside 
T are 0. Therefore the restriction of u to T' is 0. From this we deduce that the sum in 
(I4.57|l is direct, and that: 


(4.58) 

(4.59) 


= Ag(/"), for/c^n-1, (4.60) 

B”-^(/") = A”-^(/") © j;”-i(T). (4.61) 
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Recalling that A has the extension property, we see in particular that: 

dim(I^) = dim A^(r) + ^ dini£;'=(r), (4.62) 

dim T—k+1 

= 4(r)+ E (dim4'(r) +dim i?^(T)), (4.63) 

dim T^k+1 dimT—fc+1 

= ^dimRo^(T), (4.64) 

T 

where T runs through the set of faces of This proves that B has the extension 
property by Proposition^^ The sequence exactness in the form (Ell) has already been 
proved, see (|4.45|) and Lemma Td.91 Thus we have proved that R is a compatible finite 
element system. It is minimal containing A by application of Proposition 14.81 

□ 

If we want degrees of freedom for this compatible FES that yield commuting inter¬ 
polators, we can take the ones defined in [H Proposition 2.8]. 

4.4 Dimension count for “small pleasures” elements 

Neither the construction of Arnold and Awanou nor the one of Cockburn and Qiu 
resolves the question that seems most natural to us, namely, defining the minimal com¬ 
patible FES on hypercubes containing the sequences 

iP^A°(/") ^ VrA^{I^) ^ ... ^ nA”(/"). (4.65) 

The preceding general construction provides such a mcFES. Being neither the serendip¬ 
ity elements, nor the tiniest elements in the Cockburn-Qiu sense, we refer to these as 
the “Small Pleasures” elements. We provide a dimension count of the Small Pleasures 
elements with the following result. 

Lemma 4.11. 

Proof. We have an exact sequence: 

0 —^ 7^r-|-fe-l-^o(d") —^ ’Pr-|-fc-2A()(/") —^ . . . 

... ^ Vr+iAo~‘^{n VrAo~\n dVrA'^-\r) ^ 0. (4.67) 

This gives: 

k-l 

dimdVrA^-^I^) = ^(-l)Mimn+*Ag-i-*(/”). (4.68) 

i^O 

We also have: 

kerd|iP^Ag(/”) = dn+iAg-i(/”). (4.69) 

Hence, by the same formula as previously: 

k-l 

dimkerd|iP^A^'(/”) = ^^(-1)* dimlP^+i+^Ag-^-^(J”). (4.70) 

2=0 
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r = 4 

r = 5 

r = 6 

r = 7 

r = 8 

r = 9 

r = 10 

k = 

0 

0 

1 

4 

10 

20 

35 

56 

k = 

1 

11 

27 

54 

95 

153 

231 

332 

k = 

2 

45 

81 

133 

204 

297 

415 

561 

k = 

3 

35 

56 

84 

120 

165 

220 

286 


Table 1: Some computed values for dimi3g(/^) where B is a minimal compatible finite 
element system containing (14.651) . 


Therefore: 


fe-i 

dimH'=BrA5(J”) = ^(-ir(dim ^r+l+i - diiniP.+^Ag-i-^/")), (4.71) 

k-1 

= ^(-l)XdimiP,+i+,_2(„_fc+i+.)A"-'=+i+X/")- (4.72) 

dimiP,+,_2(„_fe+i+i)A"-'=+i+*(J")), (4.73) 

k-1 

= dimB,_i_,_2(n-fc)A"-'=+i+*(/"). (4.74) 

2=0 


Now, following [21 Section 3.2], the Koszul sequence on 

0 -> Uth-'- ^ Bt+iA'‘-i ^ ... ^ Bt+fc_iA"-'=+i ^ /cBt+fc_iA"-'=+i ^ 0, (4.75) 

is exact. We take t = r k — 2n, and deduce: 


k-1 

^(-l)MimBi+fc_i_.A'^-^-+i+* 
2 = 0 


dim/tBt+fe_iA” 

V k-1 j n-k )’ 

/r + 2k — n — 1\ /r + k — n — 1\ 
k-1 j[ n-k )' 


(4.76) 

(4.77) 

(4.78) 


This concludes the proof. 


□ 


Using (14.6611 and (13.811 . we can compute the dimension of Bg(/") where B is a 
minimal compatible finite element system containing (14.651) . We give some some sample 
values in Table [TJ These spaces will be investigated further in future work. 
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